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ELECTRICAL  MEASUREMENTS  OF  FERROELECTRIC 
CERAMIC  RESONATORS 


Preface 


This  study  originated  from  the  emergence  of  piezoelectric  ceramics 
as  candidates  for  high  frequency,  low-cost  resonators,  particularly  for 
telecommunications  applications.  It  was  realized  that  the  current  genera¬ 
tion  of  measurement  standards  were  inadequate  for  the  characterization 
precision  required.  Accordingly,  a  new  approach  to  the  problem  was 
developed,  and  is  detailed  in  this  report 
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Abstract 


Future  uses  of  modem  ferroelectrics  include  thin-fiSm  resonant 
microstructures  for  timing,  sensing,  and  actuation.  The  frequency  range  of 
interest  extends  from  nearly  DC  through  the  UHF  band.  At  high  frequencies 
especially,  it  is  difficult  to  obtain  precise  material  parameters  from 
broadband  electrical  measurements.  We  have  developed  highly  accurate 
lumped  equivalent  networks  for  these  frequencies,  for  the  two  canonical 
cases  of  thickness-  and  lateral-field  excitation.  Measurements  on  these 
circuits  permit  extraction  of  the  complex  dielectric,  piezoelectric,  and 
elastic  parameters  of  the  material.  In  this  report  we  discuss  both  traditional 
and  modern  electrical  network  descriptions  of  piezoelectric  resonators;  we 
also  treat  how  they  are  modified  by  the  conditions  of  excitation  and 
inclusion  of  loss  mechanisms.  The  more  general  networks  are  then 
specialized  to  yield  lumped  circuit  versions  tailored  to  the  particular 
attributes  of  ferroelectrics:  high  coupling  and  moderate  loss.  A  numerical 
example  is  worked  by  means  of  illustration. 


Introduction 


Modem  ferroelectric  materials  are  being  developed  for  a  variety  of 
mixed-effect  microdevices  that  utilize  the  elastoelectric  coupling  stemming 
from  the  piezoeffect  These  devices  will  be  required  to  be  uniform  in  their 
behavior  to  an  extent  unimagined  in  the  past;  this,  in  turn,  necessitates  the 
ability  to  characterize  and  measure  precisely  the  material  properties  of  the 
substances  from  which  they  are  fashioned.  The  properties  in  question  are 
the  dielectric,  piezoelectric,  and  elastic  constants;  the  measurements  are 
electrical  in  nature,  and  the  characterization  requires  highly  accurate 
equivalent  electrical  networks,  valid  over  a  wide  frequency  range,  to 
describe  the  physical  processes  involved.  Traditional  equivalent  circuits 
suffer  from  a  variety  of  shortcomings  that,  e.g.,  prevent  adequate  resonator 
characterization  for  the  purpose  of  extracting  material  coefficients,  as  well 
as  use  of  electronic  frequency-temperature  compensation  in  oscillator 
(TCXO)  applications.  The  network  descriptions  presented  subsequently  are 
consequences  of  the  piezoelectric  nature  of  poled  ferroelectrics.  We  begin 
with  a  brief  discussion  of  the  effect  and  its  uses.  This  is  followed  by  a 
consideration  of  traditional  equivalent  circuit  models  and  modern 
immittance  measuring  instruments.  Finally,  the  new  broadband  models  are 
introduced  and  explored,  with  a  numerical  example. 
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Piezoelectric  Resonators  and  Devices 


The  piezoelectric  effect  was  discovered  by  the  brothers  Curie  in  1880, 
and  remained  a  curiosity  until  the  early  1920s  when  its  presence  in  quartz 
was  utilized  to  realize  crystal  resonators  for  the  stabilization  of  oscillators, 
thereby  launching  the  field  of  frequency  control.1  With  the  introduction  of 
quartz  control,  timekeeping  moved  from  the  sun  and  stars  to  small,  man¬ 
made  sources  that  exceeded  astronomy-based  references  in  stability.  Since 
then,  devices  based  on  piezoelectricity  have  expanded  dramatically  in 
application.  The  quartz  resonator  has  continued  to  evolve  to  become  a 
device  capable  today  of  precision  one  million  times  greater  than  the 
original,  and  serves  as  well  as  the  “flywheel”  in  atomic  frequency 
standards,  which  make  frequency  the  most  accurate  entity  known. 

Piezoelectric  devices  are  additionally  used  as  transducers  in 
applications  from  telephone  speakers  to  sonar  arrays;  a  newer  and  rapidly 
burgeoning  area  of  utilization  is  the  integral  incorporation  of  mechanical 
actuation  and  sensing  microstructures  into  electronic  chips.  These 
microelectromechanical  structures/systems  (MEMS)  promise  signal 
sensing,  processing,  and  outputting  features  unattainable  by  exclusively 
electronic/photonic  means.3'5  In  many  traditional  applications,  applied 
voltages  conform  to  the  norms  for  electronic  circuits;  accompanying 
mechanical  displacements  are  then  usually  nanometers  or  less.  Newer 
configurations  available  today  for  actuators  often  have  much  larger 
amplitudes  at  comparable  voltages.  For  high  frequency  stabilization  and 
control,  the  configuration  of  choice  is  the  laminated  plate,  where  thin-film 
lamina  of  ferroelectric  ceramics  or  piezoelectrics  such  as  aluminum  nitride 
or  zinc  oxide  serve  as  the  resonant  membrane,  or  where  they  are  used  to 
drive  piezoelectrically  inert  layers  of  silicon.6'12  Figure  1  shows  examples  of 
modern  piezoresonator  structures. 

One  of  the  most  appealing  aspects  of  piezoelectricity  for  modern 
applications  is  the  compelling  immediacy  and  simplicity  of  the  transduction 
mechanism.  Micro-  and  nano-electronics  are  built  upon  the  behavior  of 
charged  species  subjected  to  electric  fields;  the  extreme  miniaturization  of 
these  structures  is  owed  in  large  measure  to  their  ability  to  take  on  a 
capacitor-like  form.13  Bulky  inductors  are  eliminated,  and  thin,  planar 
electrodes  introduce  the  electric  fields  to  operate  the  circuits.  These  fields 
also  provide  the  forces  required  to  drive  mechanical  motions  in  a 
piezoelectric  device.  Elastic  field/mechanical  considerations  are  therefore 
incorporated  into  the  operations  of  modern  electronic  components  in  a 
clean,  efficient,  and  very  direct  manner,  by  making  use  of  the  voltages 
resident  on  the  chips,  via  the  piezoeffect. 
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Devices  such  as  MEM  structures  provide  otherwise  unavailable 
capabilities  yet  are  embedded  in  microelectronic  circuitry.  The  full  utility  of 
these  components  cannot  be  attained  until  their  behavior  is  accurately 
characterized  in  terms  compatible  with  their  operating  environment,  viz., 
with  admittance  matrices,  and  realized  as  equivalent  electrical  networks. 
This  portion  of  the  work  is  currently  in  an  early  stage  of  development  for 
MEMS  devices  in  general.  In  the  present  paper  we  make  use  of  lossy  plate 
resonator  configurations  and  determine  the  proper  network  equivalents 
from  which  precise  materia!  constant  values  may  be  ascertained. 

St  is  interesting  to  note  that,  in  the  19th  Century,  when  mechanics 
reigned  supreme,  luminaries  like  Kelvin  and  Maxwell  “explained”  Sight 
phenomena  in  terms  of  mechanical  models  such  as  coupled 
gyroscopes.14,15  Today,  electromagnetism  and  electronics  reign,  and  the  art 
of  mechanics  seems  in  many  respects  an  all-but-forgotten  discipline.  To 
present  the  workings  of  modem  mechanical  structures  such  as 
microactuators  and  resonators  to  the  “real  world,”  we  now  make  use  of 
equivalent  electrical  circuits! 


Measuring  and  Characterizing  Equipment 

The  properties  of  ferroelectric  materials  for  electronic  applications 
are  presently  characterized  by  a  variety  of  modalities.16  These  include 
Scanning  Tunneling  Microscopy,  Atomic  Force  Microscopy,  Electron 
Microscopy  (SEM,  TEM,  STEM),  X-Ray  (Synchrotron,  XRF,  PIXE),  Surface 
(AES,  XPS,  SIMS,  RBS,  GDMS,  TDMS),  and  Spectroscopic  Techniques 
(Raman,  NMR,  ESR).  Relatively  little  attention  has  been  given  to 
determinations  of  phenomenological  characteristics  of  these  substances 
by  modern,  computer-driven  automatic  network  analyzer  and  similar,  smart 
instrumentation  metrology;  this  is  briefly  discussed  below.  These 
instruments  are  becoming  increasingly  available  for  relatively  modest 
capita!  investments,  and  represent  an  efficient  method  of  parameter 
extraction. 


Network  Analyzers 

Equivalent  circuits  of  electromechanical  systems  such  as 
piezoelectric  resonators  serve  as  a  convenient  means  of  describing 
physics  and  associated  material  parameters  in  a  form  amenable  to 
engineering  treatment  and  measurement.  Today,  computer-based 
electronic  instrumentation  exists  for  making  highly  accurate,  broadband 
measurements  of  immitance  and  scattering  parameters  of  networks  as 
function  of  frequency.  When  applied  to  the  measurement  of  ferroelectric 


resonators  the  result  is  extremely  cost-effective.  The  resonator  can  be 
quickly  and  precisely  evaluated  as  a  working  device,  and  its  performance 
judged  against  specifications  to  insure  reproducibilities  at  an 
unprecedented  level;  the  resonator  can  as  well  serve  instead  as  a  test 
structure  for  the  accurate  determination  of  the  properties  of  the  materials 
from  which  it  is  fashioned.  We  concentrate  on  the  first  aspect  in  the  paper; 
the  second  will  be  treated  in  a  sequel.  Regardless  of  the  purpose  of  the 
measurements,  the  method  depends  upon  the  availability  of  accurate 
circuit  models  characterizing  the  physics  of  the  resonating  structure. 

Once  an  accurate  equivalent  circuit  is  obtained,  a  programmable 
network  analyzer,  e.g.,  HP  8753C  (300  kHz  -  3  GHz)  or  HP  4195A  (10  Hz-  500 
MHz)  can  be  used  to  obtain  the  circuit  element  values.17  These  are  then 
related  to  the  material  values.  Alternative  measurement  methods  and 
techniques,  and  governing  standards  are  described  in  the  literature.18'24 


The  Butterworth-Van  Dyke  Circuit 

The  Butterworth-Van  Dyke  (BVD)  equivalent  circuit25,26  of  a  resonant 
system  driven  by  exposure  to  a  capacitive  electric  field  has  been  exploited 
for  many  years  to  represent  piezoelectric  resonators.27'34  The  traditional 
BVD  circuit  is  given  in  Fig.  2a.  A  distinction  must,  however,  be  made 
according  to  the  direction  of  the  exciting  field;  the  distinction  is 
represented  by  the  differences  in  Figs.  2a  and  2b.  These  will  be  considered 
in  a  later  section. 

One  of  the  best  expositors  of  the  use  of  equivalent  circuits,  and  their 
applications  to  both  microscopic  and  macroscopic  physics  and  material 
science  problems  was  Prof,  von  Hippel.27,28  He  was  also  an  exponent  of 
having  people  of  differing  technical  disciplines  able  to  translate  freely 
between  the  terminology  and  techniques  of  physics,  engineering,  and 
material  science.  We  seek  here  to  foster  this  same  facility  by  the  application 
of  electric  circuits  to  the  determination  of  material  coefficients  and 
performance  parameters  of  ferroelectric  components. 

The  BVD  circuit  of  Fig.  2a  leads  to  the  space  curve32  depicted  in  Fig. 
3,  where  the  real  (Gp)  and  imaginary  (Bp)  parts  of  the  complex  admittance 
function  are  plotted  against  frequency  variable,  f.  A  number  of  resonator 
critical  frequencies  are  indicated;  for  example,  the  resonance,  fR,  and 
antiresonance,  fA,  frequencies,  which  occur  at  the  zero  phase  points, 
where  the  susceptance,  Bp,  vanishes.  These  particular  frequencies  are  of 
considerable  importance  as  they  determine  such  things  as  bandwidth  in 
filter  applications,  and  the  normal  frequency-reactance  operating  range  of 
an  oscillator.  Projections  of  the  space  helix,  representing  the  conductance, 
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Gp,  and  susceptance  as  functions  of  frequency  are  used  to  infer  certain  of 
the  circuit  parameter  values.  Another  associated  function,  the  parallel 
capacitance,  Gp  =  Bp/co,  is  shown  in  Fig.  4  in  normalized  form,  with  Cp/C0 
plotted  against  Q  =  o»/(L1  C^.  By  measuring  the  hyperbola  constant,  and  the 
extrema,  the  BVD  elements  may  be  found.  Over  the  years  it  has  been  found 
that  the  simple  BVD  circuit  often  gives  an  inadequate  representation  when 
high  precision  measurements  are  made,  particularly  with  ceramics,  but 
even  with  quartz  resonators.20,35  For  these  instances  more  detailed 
equivalent  networks  are  invoked.  In  this  paper  the  extended  circuits  given 
stem  from  more  detailed  considerations  of  the  exact  acoustic  transmission 
line  description,  which  is  then  recast  in  lumped  parameter  form.  A  list  of 
symbols  used,  units,  and  dimensions  is  contained  in  Appendix  1. 


Simple  Thickness  Modes 

In  order  to  treat  our  subject  without  undue  complication,  we  assume 
that  the  structure  under  consideration  is  a  piezoelectric  plate  resonator 
having  an  unbounded  lateral  extent;36-41  for  all  practical  purposes  this  is 
achieved  for  the  thin-film  resonators  mentioned  earlier,  which  typically 
have  diameter/thickness  ratios  of  100  or  more.  This  assumption  permits 
use  of  a  one-dimensional  approximation,  and  leads  to  simplified 
expressions  for  the  circuit  parameters  in  terms  of  the  material  constants  of 
the  structure.  In  more  complicated  situations,42,43  the  circuit  parameters  can 
be  extracted  in  a  similar  fashion  as  well;  our  purpose  here  is  to  introduce 
improved,  yet  still  not  unduly  complicated,  networks  for  ferroelectric 
resonators. 

We  say  that  plate  resonators  execute,  in  the  one-dimensional 
approximation,  simple  thickness  modes  of  vibration.  These  are  three  in 
number,  one  is  the  thickness  extensional  mode,  the  two  others  are  the 
thickness  shear  modes;  each  mode  is  comprised  of  a  family  of  harmonics, 
representing  varying  numbers  of  phase  reversals  along  the  thickness 
direction.  In  this  paper  we  consider  one  mode  of  vibration;  the  extension  to 
three  modes  is  straightforward,  as  is  the  treatment  of  coupling  of  the 
simple  thickness  modes.  For  these  plate  modes,  it  is  found39  that  the 
capacitance  ratio,  r  =  C0/C1f  of  the  BVD  circuit  is  related  to  the  piezoelectric 
coupling  factor,  k,  by  the  relation:  r  =  7t2/8k2.  Appendices  2-4  discuss  further 
the  use  and  meaning  of  k. 
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Thickness  and  Lateral  Excitation 


A  plane  acoustic  wave  traveling  in  a  piezoelectric  medium  is 
accompanied  by  an  electric  field  parallel  to  the  direction  of  propagation. 
The  field  is  generated  by  the  wave  motion,  and  exists  for  all  directions  of 
propagation  in  the  case  of  media  of  general  anisotropy;  in  media  of  greater 
symmetry,  the  field  may  not  exist  for  certain  directions  of  travel.  In  this 
situation,  the  waves  cannot  be  transduced  by  imposed  fields  along  these 
directions,  but  it  may  be  that  they  can  be  driven  by  laterally  directed  fields, 
if  an  appropriate  piezoelectric  constant  is  present.  For  thickness  modes, 
the  direction  of  travel  is  along  the  thickness  of  the  plate.  The  usual  case  of 
transduction  is  via  a  pair  of  electrodes  placed  on  the  major  surfaces  of  the 
plate,  to  produce  a  thickness-directed  driving  field;  this  is  known  as 
thickness  excitation  (TE).  Driving  fields  impressed  from  the  side  produce 
lateral  excitation  (LE).  Figure  5  depicts  both  situations.  In  the  LE  case  the 
driving  field  is  orthogonal  to  the  wave-produced  field  in  the  thickness 
direction;  this  situation  is  appropriate  to  the  BVD  circuit  of  Fig.  2a.  Lateral 
excitation  generally  permits  greater  design  latitude  because  the  azimuth 
angle  provides  an  additional  degree  of  freedom.  In  ceramics,  the  thickness 
direction  often  coincides  with  the  poling  direction,  and  the  lateral  direction 
is,  therefore,  in  the  isotropic  basal  plane  where  the  azimuthal  angle  is  of  no 
consequence.  The  TE  situation,  where  the  driving  and  wave-produced  fields 
are  parallel,  is  shown  by  the  theory31,34,39,44^6  to  be  characterized  by  the 
BVD  circuit  of  Fig.  2b;  the  presence  of  the  negative  capacitor  is  a 
manifestation  of  the  interaction  of  the  driving  and  wave-produced  fields. 
Thickness  excitation  is  the  overwhelmingly  predominant  form  of  excitation 
for  piezoelectric  plates. 

In  Fig.  6  the  BVD  circuits  for  LE  and  TE  are  contrasted.  It  is  seen  from 
the  figure  that  the  capacitance  ratios  r  and  r’  differ  by  unity.  A  number  of 
important  electrical  attributes,  e.g.,  the  frequency  difference  (fA  -  fR),  of  the 
networks,  are  related  to  the  capacitance  ratio.  Depending  upon  the  size  of 
the  piezocoupling  factor,  k,  the  use  of  the  traditional  BVD  network  of  Fig.  2a 
(pertinent  to  the  LE  case),  when  the  circuit  of  Fig.  2b  is  called  for  (in  the 
usual  case  of  TE),  can  result  in  a  significant  error.  This  will  be  seen  in 
greater  detail  in  a  later  section  when  a  ferroelectric  ceramic  is  used  as  an 
example.  For  now,  consider  the  case  of  a  poled  ceramic  where  k  =  51%. 
From  the  relation  given  above,  r  =  rc2/8k2,  we  find  that  r  =  4.74,  and  r’  =  3.74; 
the  difference  is  greater  than  20%,  and  is  therefore  quite  significant  For  a 
low  coupling  substance  such  as  AT-cut  quartz,  where  k  =  8.80%,  r  and  r’ 
differ  by  only  0.6%;  this  is  frequently  of  no  consequence. 
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In  Fig.  7  is  plotted  the  reactance  of  a  resonator  of  a  high-coupling 
material  versus  frequency  for  both  TE  and  LE  methods  of  excitation.  The 
relative  shift  in  resonance  and  antiresonance  frequencies  in  the  two  cases 
is  readily  apparent.  Resonance  for  LE  coincides  with  the  TE  antiresonance. 
The  effect  diminishes  with  the  square  of  the  harmonic,  as  does  the  pole- 
zero  separations;  this  is  a  consequence  of  the  diminution  of  effective 
piezocoupling  with  overtone. 


Exact  Lossless  Network 

Synthesis  of  the  network  exactly  realizing  Yin,  the  input  admittance 
of  a  plate  resonator  executing  simple  thickness  modes,  discloses  that  the 
individual  circuit  elements  can  be  given  simple  interpretations  in  terms  of 
the  physical  processes  within  the  mechanically  vibrating  structure.  This 
leads  to  analog  circuit  descriptions  of  piezoresonators.  The  complete 
analog  network  for  the  TE  case  is  shown  in  Fig.  8.  The  transmission  line 
represents  propagation  of  acoustic  waves  within  the  piezoelement. 
Extension  to  the  case  of  the  composite  resonator48,49  is  shown  in  Fig  9.  The 
network  representation  follows  by  inspection  because  of  the  modular 
nature  of  the  circuit  parts.  Composite  resonators  are  important  for 
determining  the  material  coefficients  of  nonpiezoelectric  materials,  and  as 
substrates  for  thin-film  resonators  (Fig. 1b).  It  is  seen  In  Fig.  8  that  the 
circuit  accommodates  electrode  mass  deposited  on  the  major  surfaces  of 
the  plate.  In  the  figure  the  electrode  is  considered  to  be  thin  enough  that 
acoustic  wave  propagation  effects  therein  can  be  neglected,  and  the  effect 
of  mass  loading  wholly  ascribed  to  inertia,  which  is  represented  by  an 

50-53 

inductor  at  each  mechanical  port. 

If  the  electrode  mass  loadings  on  each  surface  of  the  plate  are  equal, 
the  circuit  structure  of  Fig.  8  is  completely  symmetric,  and  it  may  be 
bisected.  The  result  is  seen  in  Fig.  10.  A  switch  has  been  added  as  well,  to 
depict  the  effect  on  the  circuit  of  changing  from  TE  to  LE.  Exact 
expressions  for  the  input  admittance,  Yin,  of  a  single  excited  mode  in  the 
one-dimensional  approximation,  when  driven  in  TE  and  LE,  are: 

Yin(TE)  =  j  co  C0  /  [1  -  k2T(X)3  ;  Yin(LE)  =jcoC0  [1  +  k2T(X)3 

where  T(X)  =  tan(X)/(X),  and  X  =  Kh.  The  parameter  k  =  co /v  is  the  acoustic 
wavenumber.  Piezoelectric  transformer  turns  ratio,  n,  appearing  in  Fig.  10, 
is  found  from  the  relation  n2  =  (eA/2h)  *  2  C0k 2f1/Y0;  see  Appendix  1. 
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Inclusion  of  Loss 


54-64 


The  exact  circuit  of  Fig.  8  is  for  the  lossless  situation.  We  may 
include  the  phenomenological  effects  of  loss  by  two  simple  expedients. 
Internal  friction  within  the  piezo  material  is  accounted  for  by  making  the 
elastic  stiffness  complex:  cE  -»  c*  =  cE  +  jcori;  the  imaginary  part  is  the 
acoustic  viscosity  term.58  This  has  the  effect  of  making  the  velocity  and 
wavenumber  complex,  and  leads,  in  the  lumped  BVD  description,  to  the 
presence  of  the  R1  term.  The  second  provision  is  for  dielectric  and  ohmic 
losses.  We  make  the  dielectric  permittivity  complex:  ss  ->  s*  =  s’  -  js”,2and 
combine  ohmic  and  dielectric  loss  terms;  the  total  effect  is  additive,  so 
that  the  total  effective  a  is  the  ohmic  conductivity  a  plus  the  dielectric 
losses,  cos”.  In  the  traditional  BVD  circuit  these  losses  appear  as  a  shunt 
resistor  across  the  static  capacitor  C0;  (this  can  be  thought  of  as  a  lossy 
capacitor  C*).  When  the  LE  version  of  the  BVD  circuit  is  taken,  the  further 
addition  of  the  shunt  resistor  constitutes  the  ceramic  resonator  equivalent 
circuit  almost  universally  used  at  low  frequencies.35  Viscous  and  dielectric 
mechanisms  are  separately  depicted  in  Fig.  11.  The  motional  time  constant 
is  tt  =  R-,0.,  =  t|/c,  where  c  is  the  lossless  piezoelectrically  stiffened  elastic 
stiffness  for  the  mode  in  question  (Appendix  4).  The  “static”  time  constant 
is  t0  =  R0C0  =  s’/  (a  +  cos”),  which  reduces  to  x0  -  s’/  a  in  the  DC  limit.  Both  x-, 
and  x0  are  real  numbers. 

In  general,  both  s’  and  s”  are  functions  of  frequency,  while  s’(co)  and 
s”(oo)  are  themselves  interrelated  by  the  Cauchy-Riemann  equations.  In 
such  a  case,  the  network  for  C0*  is  a  parallel  connection  of  a  real  C0  and  a 
real  G0,  both  of  which  are  functions  of  frequency.  Alternatively,  depending 
upon  the  frequency  variation  (Cole-Cole  relaxation,  resonance,  etc.),  the 
circuit  might  be  given  a  more  complicated  topological  form,  but  with 
elements  having  frequency-invariant  values.  A  variation  of  s*  with 
frequency  also  adds  to  the  dispersion  in  the  acoustic  transmission  line  that 
arises  from  the  lossy  stiffened  elastic  stiffness.  Additional  details  are 
contained  in  Ref.  [65]. 

Lumped  Circuits 

Passage  from  the  distributed  transmission  line  description  to  lumped 
parameters  is  accomplished  by  using  the  partial  fractions  expansion  of  the 
tangent  function  describing  the  transmission  line,  as  shown  in  Fig.  12. 
Loss  is  included  by  making  the  argument  of  the  tangent  function  complex, 
as  described  earlier.  The  individual  partial  fractions  then  each  have  lumped 
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network  equivalents  as  series  RLC  circuits,  and  represent  the  harmonic 
plate  resonances. 

The  bisected,  exact,  transmission  line  circuit  model  characterizing  a 
lossy  plate  resonator  is  shown  in  Fig.  13.  Included  here  is  the  extension  to 
the  case  where  the  electrode  is  of  a  thickness  where  its  wave  transmission 
time  cannot  be  neglected,  and  the  electrode  must  itself  be  represented  by  a 
transmission  line.  The  figure  shows  the  presence  of  ohmic  and  dielectric 
conductivities  modeled  by  the  shunt  conductance  G0  =  1/R0,  as  in  the  BVD 
case.  Of  particular  note  is  the  presence  of  a  negative  G0  shunting  the 
negative  C0  of  the  TE  circuit;  the  negative  G0  and  C0  arise  in  the  lossy  case 
for  the  same  reason  as  the  negative  G0  in  the  lossless  case:  the  interaction 
of  the  driving  and  wave-produced  fields. 

When  one  simplifies  Fig.  13  appropriately,  one  obtains  the  network  of 
Fig.  14,  the  new  high  frequency  TE  resonator  equivalent  circuit  pertinent  to 
high  coupling  ceramic,  and  other  ferroelectric  materials.  For  LE,  one  simply 
short-circuits  the  negative  elements  to  arrive  at  the  appropriate  network 
topology.  As  regards  the  numerical  circuit  values  pertinent  to  the  two 
excitation  types,  we  note  that  values  of  k  for  TE  and  LE  are  mutually 
independent,  stemming  from  different  piezoelectric  coefficients;  generally, 
the  electrode  geometry  is  also  different,  and  therefore  the  two  types  of 
piezoelectric  driving  are  represented  not  only  by  different  circuit 
topologies,  but  also  by  different  numerical  values.  The  simplifications 
necessary  to  yield  Fig.  14  from  the  distributed  analog  network  are  the 
following: 


1)  a  single  resonant  harmonic  is  retained,  and  modeled  by 
RM,  Lm,  Gm;  M  is  an  odd  number  representing  the  harmonic. 

2)  the  effect  of  all  other  harmonics  is  subsumed  into  another 
capacitor,  CR,  shunting  the  RM,  LM,  CM  branch. 

3)  electrode  inertia  is  modeled  by  an  inductor  Le. 

The  more  accurate  circuit  of  Fig.  14,  or  the  simpler  version  of  Fig.  2c 
for  the  situation  when  CR  can  be  neglected,  is  to  be  contrasted  with  the 
network66'69  usually  given  for  ceramic  resonators;  viz.,  an  LE  BVD  with 
shunt  G0.  If  necessary  to  expand  the  frequency  measurement  regime  to 
more  than  one  harmonic,  the  circuit  of  Fig.  14  may  be  augmented  by  the 
placement  of  additional  series  RLC  arms  shunting  that  shown,  with  a 
concomitant  decrease  in  the  value  of  CR. 
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Circuit  Element  Values 


We  assume  a  thickness-excited  (TE)  piezo  resonator  of  thickness  2h, 
coated  with  coextensive  electrodes  of  area  A,  mass  density  pe,  and 
thickness  he  on  the  major  plate  surfaces.  The  plate  is  of  mass  density  p; 
complex  dielectric  permittivity  and  total  effective  conductivity  in  the 
thickness  direction  are  s  and  a,  effective  piezoelectric  constant  is  e,  and 
the  stiffened  elastic  constant  and  viscosity  governing  the  modal  velocity 
and  attenuation  are  denoted  by  c  and  r|.  The  requisite  values  of  s* ,  a,  e,  t| 
and  c  are  obtained  from  the  material  tensor  schemes  by  straightforward 
methods.39,44 

In  terms  of  these  quantities,  and  the  definitions  in  Appendix  1,  the 
circuit  elements  of  Fig.  14  are,  to  our  level  of  approximation: 


C0  =  e5  A/2h  R0  =  1/G0  =  xo/C0 

cM  =  8  k2  c0  /mV  cr  =  c, (71^/8  -  i/m2) 


Lm  =  h2/2v2C0  k2 


Le  =  ph2/v2C0  k2  =  2 p  Lm 


k2  =  e2/s’c 


More  generally,  k2  is  complex,  which  renders  the  circuit  elements  complex. 
The  effects  can  virtually  all  be  subsumed  into  the  same  circuit  form  with 
modified  values;  e.g.,  LM  and  CM  each  contain  small  series  resistances, 
which  would  be  added  to  RM. 


The  curious  relation  Le  =  2|i  LM  arises  for  the  following  reason: 
mechanical  displacement  along  the  thickness  direction  is  sinusoidal,  with 
antinodes  at  the  plate  surfaces;  an  electrode  mass  lumped  at  the  surface 
has  twice  the  influence  of  equivalent  mass  distributed  throughout  the 
thickness  because  the  latter  must  be  weighted  by  the  square  of  the 
distribution  function. 


Ferroelectric  Materials.  Properties,  and  Preparation 

References  [70]-[85]  give  data  on  the  measured  constants  of  many 
ferroelectric  materials,  as  well  as  various  properties  and  methods  of 
preparation.  In  Fig.  15  is  shown  a  qualitative  depiction  of  the  loci  occupied 
by  various  classes  of  piezoelectric  materials  in  the  quality  factor  vs 
capacitance  ratio  plane.  Ceramics  occupy  a  unique  and,  from  a 
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measurement  standpoint,  a  somewhat  troublesome  position,  possessing 
very  small  r  values  (high  coupling  factors)  and  low  to  moderate  Q  values. 
Because  of  this  position,  loss  considerations  are  very  important,  as  is  the 
distinction  between  TE  and  LE  methods  of  excitation,  and,  consequently, 
use  of  the  appropriate  equivalent  network. 

Tables  S,  II,  and  III  provide  the  matrices  of  the  permittivity, 
piezoelectric  coefficient,  and  elastic  stiffness  for  materials  in  crystal  class 
6mm,  which  includes  the  piezoelectric  ceramics,  discussed  in  Appendix  5. 
Values  of  ct  and  r|  for  ferroelectrics  in  general,  and  piezoelectric  ceramics  in 
particular,  are  quite  scarce  in  the  literature.  This  may  be  at  least  partially 
attributable  to  high  variability  from  the  lack  of  sufficient  reproducibility  in 
manufacture;  the  existence  of  wide  variations  in  loss  with  temperature,  and 
hysteretic  effects  may  also  be  responsible  for  the  dearth  of  these  data.  Part, 
at  least,  of  the  reason  is  due  to  the  use  of  improper  or  inaccurate  equivalent 
circuit  representations  by  which  to  characterize  these  substances.  This 
results  in  circuit  element  values  that  are  functions  of  frequency,  and  which 
cannot  be  related  in  an  obvious  and  reproducible  manner  to  material 
parameters. 


Numerical  Example 

We  give  as  an  example  the  circuit  parameters  and  electrical  input 
behavior  of  a  thin-fiim,  thickness-mode  ceramic  resonator  as  might  be 
fabricated  by  sol-gel,76’85  MOCVD,  or  laser  ablation  methods  for  operation  in 
the  VHF-UHF  frequency  bands.  The  material  chosen  is  a  modified  lead 
titanate.83  Piezoelectrically  stiffened,  lossless  acoustic  velocity  is  v  =  4400 
m/s;  piezocoupling  is  k  =  51%;  and  dielectric  constant  in  the  thickness 
direction  is  eT33/s0  =  263.  We  choose  specifically  a  fundamental  mechanical 
response  frequency  of  f,  =  250  MHz;  the  resonator  thickness  required  is  2h 
=  v/2f,  =  8.8  urn.  Mass  density  is  taken  as  7.52  IQ3  kg/m3.  To  obtain  the 
permittivity  ss,  necessary  to  calculate  C0,  one  would  compute  es33  =  sT33  - 
e33d33  -  2e31d31.  This  follows  from  the  relation  As  =  ed’  in  Appendix  3  applied 
to  substances  with  transverse  isotropy,  using  the  entries  in  Table  II.  Since 
values  of  e31  and  e33  were  not  published  in  Ref.[82J,  we  make  the  reasonable 
approximation  that  ss  »  sT33  (1-k2).  This  yields  ss  *  184.59  s0  «  0.7399  eT33  « 
1.723  nF/m.  The  piezo-transformer  turns  ratio  squared  is  therefore  n2  « 
0.842.  We  also  assume  a  viscous  (or  motional)  time  constant^.,  of  value  T.j  = 
10'12  s  (loss  one  hundred  times  greater  than  quartz);  this  is  presently 
unattainable  at  VHF/UHF,  but  is  anticipated  with  future  advances  in  material 
science.  A  value  of  total  a  (conduction  plus  dielectric  loss)  of  10‘7  S/m  is 
used  in  the  calculations,  corresponding  to  a  static  time  constant,  x0  »  ss/ar, 
of  17.23  ms. 
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We  assume  also  a  square  electrode  patch  on  each  surface  of  the  thin 
film  with  dimensions  1  mm  by  1  mm  by  250  nm.  The  electrode  area  is  thus 
A  =  10'6  m2,  and  the  width-to-plate-thickness  ratio  is  1  mm/8.8^m  «  114,  so 
the  one-dimensional  approximation  is  valid.  The  acoustic  (mechanical) 
impedance  of  the  ceramic  for  this  particular  mode  is  Z0  =  A  pv  »  33.09 
mechanical  ohms  (kg/s). 

Assuming  further  that  the  electrode  is  of  aluminum,  with  mass 
density,  pe  of  2.70  103  kg/m3,  the  normalized  mass  loading  p  =  (electrode 
mass  per  unit  area/resonator  mass  per  unit  area)  =  pehjph  = 
(2.70x250nm)/(7.52x4.4pm)  *  2.04%.  Because  the  acoustic  velocity  in  the 
electrode  is  roughly  of  the  same  magnitude  as  that  in  the  ceramic,  this  is 
also  very  nearly  the  ratio  of  the  acoustic  lengths  of  the  electrode  and  plate 
in  the  thickness  direction.  Because  p  is  small,  it  is  therefore  permissible  to 
substitute  a  lumped  inductor  Le  for  the  transmission  line  representing  the 
electrode;  see  Fig.  13.  Traditional  treatments  of  mass  loading50'53  show  that 
the  normalized  mass  loading  approximates  closely  the  ensuing  frequency 
lowering  due  to  the  electrode  mass;  more  accurate  results  are  reported  in 
Ref.  [39].  In  our  case,  p  =  Af/f,  «  2.04%,  or,  since  f,  =  250  MHz,  the  mass 
loaded  frequency  of  the  purely  mechanical  resonance  is  reduced  to  244.9 
MHz;  we  shall  see  that  this  is  indeed  the  effect  of  the  inductor  Le, 
representing  the  mass  loading.  For  the  purposes  of  comparison,  equal¬ 
thickness  electrodes  of  chromium  (pe  =  7.19  IQ3  kg/m3;  p  *  5.45%),  and  of 
gold  (pe  =  19.32  103  kg/m3;  p  «  14.61%),  will  also  be  simulated. 

The  above  data  permit  us  to  calculate,  using  the  formulas  given 
earlier,  the  parameters  associated  with  the  equivalent  network  of  Fig.  14; 
these  are  given,  to  an  accuracy  consistent  with  our  other  approximations, 
as  follows: 

C0  *  1 95.8  pF;  k2  C0  ~  50.93  pF 
R0  =  1/G0  =  8.8  107  ohm 
Le  «  0.401  nH  (aluminum) 

«  1.07  nH  (chromium) 

*  2.87  nH  (gold) 

For  M  =  1  (fundamental  harmonic), 

Cm  =  C1  -  41.3  pF 
CR  »  9.65  pF;  CM  +  CR  «  50.95  pF 
RM  =  R1  ^  1/41.3  ohm 
LM  =  L,  »  9.82  nH 
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Figure  16a  shows  broadband  simulations  of  the  input  admittance  |Yin| 
as  functions  of  frequency,  using  the  new,  advanced  equivalent  circuit  for 
electrodes  of  gold,  chromium,  and  aluminum,  showing  the  effects  of 
differing  values  of  the  mass  loading  parameter,  p;  the  curves  for  Gp  are  very 
similar  to  those  of  |YJ.  Fig.  16b  shows  the  input  capacitance  Cp  versus 
frequency,  for  aluminum  electrodes;  the  corresponding  extrema  are  given 
in  Table  IV.  These  values  are  to  be  compared  with  those  in  Table  V, 
computed  for  the  situation  where  the  electrode  is  of  negligible  mass.  It  is 
seen  that  the  primary  effect  is  a  diminution  of  the  critical  frequencies. 

In  Table  VI  are  given  the  resonance  and  antiresonance  frequencies 
(defined  as  the  zeros  of  the  Cp  function),  for  a  number  of  electrode 
substances  of  differing  mass  loading  (p)  values.  The  resonance  frequencies 
correspond  closely  with  those  of  the  admittance  maxima  of  Fig.  16a.  The 
last  column  of  Table  V!  is  instructive,  as  well  as  of  considerable  practical 
importance.  The  frequency  difference  (fR  -  fA),  is  associated  with  the 
bandwidth  of  filters,  and  with  the  inductive  (operating)  region  of  resonators. 
It  is  determined  jointly  by  three  factors:  piezocoupling,  loss,  and  mass 
loading.  For  a  lossless  resonator  having  our  assumed  material  and 
geometrical  parameters,  without  electrode  mass  loading,  the  fractional 
frequency  difference  (fA  -  fR)/f,  is  [V(1+1/r)  -  1]  »  10.04%.  Addition  of  loss 
increases  this  to  *  11.78%;  this  is  found  from  the  first  entry  in  the  last 
column  of  Table  VS.  Mass  loading  decreases  the  fractional  difference 
slightly,  by  only  «  2.81%  per  percent  p.  To  a  good  approximation,  the  fA(i.e., 
purely  mechanical)  frequencies  may  be  obtained  from  the  relation  fA  =  (1  -p) 
fi- 


The  foregoing  has  considered  the  case  of  thickness  excitation.  To 
examine  lateral  field  excitation,  we  assume  that  the  lateral  coupling 
coefficient  is  also  51%.  In  the  interest  of  direct  comparison,  we  assume  that 
2.04%  mass  loading  is  present  as  well,  although  LE  units  typically  are 
devoid  of  coatings  on  their  major  faces.  Dielectric  coatings  are  sometimes 
used  for  frequency  adjustments,  but  will  in  the  future  be  discouraged 
because  of  aaina  effects,  and  other  disadvantages.  Table  VII  displays  the 


results  for  LE;  they  are  obtained  by  short-circuiting  the  -G0  and  -G0 
elements  in  the  TE  network.  Comparison  with  Table  IV  discloses  that  the 
responses  for  LE  are  reduced  in  magnitude;  this  follows  from  the  relations 
in  Fig.  6;  because  the  LE  capacitance  ratio  is  larger  than  for  TE,  the 
effective  coupling  is  diminished,  and  the  responses  made  weaker.  The 


Table  VIII  gives  the  results  when  both  CR  and  Le  are  set  equal  to  zero; 
this  corresponds  to  zero  mass  loading  plus  neglect  of  all  harmonics  except 
the  fundamental.  The  Table  IX  results  retain  the  finite  aluminum  mass 
loading,  but  again  have  CR=  0.  Table  X  through  Table  XVII  list  the  results  for 
eight  different  parameter  conditions  (Gp  =  maximum,  etc.),  when  the  five 
circuit  variations  discussed  above  (CR=  0,  etc.)  are  made. 

Table  XVIII  through  Table  XX  give  the  resonance,  antiresonance,  and 
fractional  difference  frequencies  for  aluminum,  chromium,  and  gold 
electrodes,  respectively.  For  those  two  circuit  conditions  having  Le  =  0,  the 
mass  loading  vanishes,  and  the  corresponding  frequencies  in  all  three 
tables  are  seen  to  agree,  as  they  must.  Table  XXI  provides  the  susceptance 
and  capacitance  slopes  at  the  inflection  points  of  the  corresponding 
functions.  Table  XXII  shows  the  frequencies  where  the  admittance  is 
maximum  for  the  complete  circuit  of  Fig.  14,  and  for  four  circuit  fragments. 


Behavior  at  Low  Frequencies 

Very  often  measurements  are  made  at  low  frequencies,  typically  1 
kHz,  in  order  to  infer  values  for  the  shunt  resistor  R0,  and  the  dissipation 
factor  Df.  This  latter  factor  is  arbitrarily  defined  as  the  quotient  (Gp/Bp)  at  1 
kHz,  and  is  often  set  equal  to  the  reciprocal  of  a  quality  factor,  Q’.  While 
qualitative  inferences  result,  the  notion  of  using  such  measurements  is  at 
best  ambiguous,  and  must  be  treated  with  some  care.  The  resonator  quality 
factor,  Q,  is  best  regarded  instead  as  the  quantity  Q  =  V(LM/CM)/RM;  this 
depends  solely  on  the  elements  of  the  motional  arm,  is  independent  of 
frequency,  and  helps  determine  the  vibrator  characteristics  of  the  unit  near 
the  resonance  frequency.  As  an  example  of  the  difficulties  that  may 
otherwise  arise,  consider  the  result  of  using  measurements  at  1  kHz  to  fix  a 
value  of  R0.  At  this  frequency  the  input  capacitance  is  Cp«  264.67  pF,  and 
the  Q’  w  166.9.  Therefore,  R0  =  Q7(2jtfCp)  «  100.4  MQ,  not  88  MQ,  the 
constant  value  approached  in  the  DC  limit.  Defining  Q’  as  (Bp/Gp)  results  in 
a  frequency  dependent  value  which  approaches  zero  (and  Df  becomes 
unbounded)  as  f->  0,  a  situation  devoid  of  practical  significance. 

Table  XXIII  gives  the  values  of  Q’  =  (Bp/Gp)  and  Df  =  1/Q’  evaluated  at 
1kHz  for  the  complete  circuit  of  Fig.  14,  and  for  four  circuit  fragments. 

In  the  DC  limit,  moreover,  one  must  distinguish  carefully  between  the 
lossless  and  lossy  cases  when  using  the  TE  version  of  the  circuit.  In  the 
lossless  case,  the  presence  of  the  -C0  element  results  in  an  all-capacitor 
network  having  an  input  admittance  of  jcoC^I  -  k2)  as  co  -*  0;  in  the  lossy 
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an  exact  model  for  a  single,  simple  thickness  mode,  and  could  be 
augmented  to  include  additional  modes.  This  would  have  added  to  the 
complexity  without  adding  to  insight.  It  could  have  also  been  modified  to 
incorporate  nonuniform  lateral  distributions  of  energy,  and  “energy 
trapping”  considerations.39  In  these  cases,  the  new  circuit  of  Fig.  14  would 
retain  its  lumped  topological  form,  but  the  element  values  would  be 
modified.  For  example,  the  presence  of  other  harmonics  of  a  single  mode 
type  is  approximately  accommodated  in  the  improved  equivalent  circuit  by 
the  provision  of  CR;  addition  of  other  modes  would  simply  increase  its 
value  in  a  calculable  manner. 

Having  extracted  the  advanced  network  of  Fig.  14,  the  question 
naturally  arises  whether  one  can  further  reduce  it  to  the  form  of  an 
elementary  BVD  circuit  (Fig.  2a)  with  shunt  resistance,  having  frequency- 
independent  elements.  This  is  not  possible,  as  network  analytical  methods 
will  show.  One  may,  however,  construct  a  rudimentary  BVD  circuit  that 
“best”  describes  the  complete  resonator  in  some  fashion.  There  are 
obviously  many  ways  to  accomplish  this;  we  give  below  one  such,  for  the 
purpose  of  comparison  with  the  more  advanced  network  of  Fig.  14. 

We  arrive  at  an  equivalent  rudimentary  (five-element)  BVD  by  making 
the  following  assumptions: 

®  G0  equals  that  of  the  advanced  circuit  as  f-»  oo 
©  C0  +  C1  equals  that  of  the  advanced  circuit  as  f 0 
These  two  criteria  fix  C0,  C1s  and  r. 
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•  Gp  equals  G0  of  the  advanced  circuit  as  f-»  0 

This  fixes  R0  =  1/G0. 

•  (^pmax  “  Cpmin)  =  Q  C1  =  F  C0  equals  that  of  the  advanced  circuit 

This  fixes  Q,  F,  E,  and  the  f,  R.,  product. 

•  fi  equals  the  arithmetic  average  of  the  following  four  values: 

f-11  =  value  at  which  the  f,  R1  product  is  satisfied  when  R.,  is 
taken  as  1/Gpmax  of  advanced  circuit 
fi2  =  (f  at  Cpmax  -fat  Cpmjn)/2  of  advanced  circuit 
f13  =  f  at  Gpmax  of  advanced  circuit 
f14  =  f  (between  Cpmax  and  Cpmin)  at  which  Cp  of  advanced 
circuit  equals  C0. 

This  fixes  the  values  of  R1  and  L,. 

These  four  frequency  values  would  be  identical  if  the  rudimentary  BVD  was 
an  adequate  equivalent  circuit  The  above  recipe  yields  the  following 
parameter  values: 

R0=  88  MQ;  C0=  195.8  pF; 

L-i  =  1 0.70  nH;  R1  =  (1/46.38)  Q;  C,  =  (CR+  CM)  =  50.95  pF 

A  comparison  of  some  pertinent  quantities  of  the  advanced  and 
simple  lumped  networks  is  given  in  Table  XXIV.  The  slope  of  the  normalized 
capacitance  curve  at  the  inflection  point  is  equal  to  -2FQ  =  -2E,  for  the 
simple  BVD  circuit.86,87  The  table  indicates  discrepancies  from  5%  to  greater 
than  20%.  Simulations  of  input  conductance  and  capacitance  versus 
frequency  for  both  the  advanced  and  simple  networks  are  shown  in  Fig.  17a 
and  b,  respectively,  in  the  vicinity  of  the  resonance.  It  is  seen  from  Fig.  17 
that  the  simple  circuit  produces  curves  that  are  qualitatively  similar  to 
those  generated  by  the  advanced  network,  but  that  the  curves  are 
sufficiently  divergent  not  to  be  adequate  for  characterizing  devices  to  the 
extent  that  will  be  required  by  future  manufacturing  specifications. 

A  comparison  of  the  curves  generated  by  the  advanced  network  of 
Fig.  14  with  those  of  the  exact  transmission  line  circuit  of  Fig.  13,  shows,  by 
comparison,  deviations  at  the  percentage  level,  and  only  in  the  immediate 
vicinity  of  fR,  the  resonance  frequency.  Some  of  the  divergence  is  due  to 
errors  in  rounding  the  element  values,  and  some  to  using  only  the  real 
parts  of  the  lumped  circuit  elements;  much  of  the  remainder  is  thought  to 
be  due  to  treating  all  the  nonresonant  harmonics  on  the  same  footing,  as 
lumped  capacitors,  instead  of  RC  combinations.  The  advanced  network, 
without  additional  modifications,  is  appropriate  and  adequate  for 
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APPENDICES 


Appendix  1.  Units  and  Dimensions 


Field  tensors: 

Mechanical  stress,  S  pascal,  Pa  =  newton/meter2  =  N/m2 

Mechanical  strain,  T  meter/meter  =  m/m  =  dimensionless 


Electric  intensity,  E  volt/meter  =  V/m 

Displacement,  D  coulomb/meter2  =  C/m2 
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Matter  tensors: 

Mass  density,  p,  pe 

kg/m3 

Permittivity,  sT,  ss 

farad/meter  =  F/m  =  C/(m-V) 

Impermeability,  pT,ps 

meter/farad  =  m/F  =  (m-V)/C 

Conductivity,  a 

siemens/meter  =  S/m 

Piezo  constants  e 

C/m2  =  N/(m-V) 

d 

mN  =  C/N 

h 

V/m  =  N/C 

g 

m2/C  =  (m-V)/N 

Elastic  stiffness,  c  ,  c 

Pa  =  N/m2  =  J/m3 

Compliance,  sE,  sD 

(Pa)'1  =  m2/N  =  m3/J 

Viscosity,  r| 

pascal-second  =  Pa-s 

Network  and  resonator  parameters  (subscript  "e"  =  electrode): 

Acoustic  impedance 

kg/s 

Z0  =  Apv  -  1/Y0 

Ze  =  Apeve  =  1/Y  e 
Acoustic  velocity 

m/s 

v  =  V(c/p) 
ve  =  V(  cjpe) 

Acoustic  wavenumber 

-i 

m 

k  -  co/v 

Ke  =  C0/Ve 

Capacitance,  C0,  Cl5 

farad 

CM,  CR 

Capacitance  ratio 

dimensionless 

r  ™  C0/C1 

Dissipation  factor 

dimensionless 

Df  =  Gp/Bp, 
at  1  kHz 

Elastic  stiffness 

pascal 

c,  ce 

Electrode  area 

meter2 

A 

Electrode  thickness 

meter 

(each) 

he 

Figure  of  excellence 

dimensionless 

E  =  Q2/r 

Figure  of  merit 

dimensionless 

F  =  Q k 
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Frequency  variables 

(0,  f 

hertz  =  Hz  =  s1 

Harmonic  number,  M 

odd  integer 

Inductance,  L,,  LM,  Le 

henry 

Mass  loading 

parameter,  p 

dimensionless 

Normalized  frequency 

Q  =  to/©!  =  f/fj 

dimensionless 

Piezocoupling  factor 
k  =  e/V(  sc) 

dimensionless 

Piezo-transformer 
turns  ratio 
n  =  Ae/(2h) 
n2  =  2C0  k^Zo 

N  N  =  C/m 

Plate  thickness 

2h 

meter 

Quality  factor 

Q  =  [^(L,/C1)]/R1 

Q’  =  1/D, 

dimensionless 

Resistance,  R0,  R^  RM 

R0  =  1/G0 

ohm 

Series  resonance 
frequency 

fj  = 

hertz 

Series  resonance 
pulsatance 
©.,  =  i/V^c^ 

hertz 

Surge  impedance 

Zi  =  ^(L,ICJ 

=  QR1 

ohm 

Time  constant 

second  =  s 

t-i  =  R1C1  =  r\lc 

(motional,  dynamic) 

To  =  R0c0  = 

e’/(< r  +  ©s”)  *  8 S/cr 

(“static”) 

Appendix  2.  Material  Constitutive  Equations 

Piezoelectricity  couples  the  mechanical  stress  (T)  and  strain  (S) 
fields  with  the  electric  intensity  (E)  and  displacement  (D)  fields.  Depending 
on  variable  choice,  four  sets  of  material  constitutive  relations  are  defined 
for  the  linear  case  of  insulating  dielectric  materials;  these  are,  in 
compressed  matrix  form: 
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are  tensors  of  rank  two  which  have  been  converted  to  six-dimensional 
vectors,  [  cE  ],  [  cD  ]  and  [  sE  ],  [  sD  ]  are  (6  X  6)  symmetric  matrices 
representing  the  elastic  stiffnesses  and  compliances  at  constant  E  and  D; 
[  £S  ],  [  eT  ]  and  [  (3s  ],  [  pT  ]  are  symmetric  (3  X  3)  dielectric  permittivity  and 
impermeability  matrices  at  constant  S  and  T,  while  [  e  ],  [  h  j,  [  d  ],  and  [  g  ] 
are  (3  X  6)  piezoelectric  matrices.  A  prime  denotes  transpose.  Alternative 
formulations,  using  the  electric  polarization  vector  [  P  ]  in  place  of  [  D  ],  are 
sometimes  used,  particularly  in  electrooptic  applications;  piezoelectric 
polarization  [  a  ]  and  [  b  ]  matrices  are  then  used.  The  matrix  brackets  [  J 
are  omitted  in  Appendix  3  and  4  for  simplicity. 

The  crystallographic  point-group  symmetry  elements  completely 
determine  the  form  and  symmetries  of  the  elastic,  piezoelectric,  and 
dielectric  matrices  that  appear  in  the  constitutive  relations.  The  specific 
values  of  the  individual  components  result,  of  course,  from  the  particular 
substance  considered. 


Appendix  3.  Matter  Tensor  Identities 
The  material  constants  are  interrelated  as  follows: 

s*  p+  =  c*  s”  -  S  ,  where  +  is  T  or  S,  *  is  E  or  D,  and  1  =  unit  matrix,  and 
e  =  ss  h  =  d  cE;  d  -  et  g  -  e  sE;  h  =  ps  e  =  g  cD;  g  =  pT  d  =  h  sD 


In  addition  to  these  multiplicative  identities,  difference  relations  exist 
between  the  dielectric  quantities  at  constant  stress  and  at  constant  strain 
and  between  the  elastic  constants  at  constant  electric  field  and  at  constant 


(cD  -  cE)  =  h’  e  =  e’  ps  e  =  h'  es  h 
(sE  -sD)  =  g'd  =  d’pTd  =  9’ sT g 
(eT  -  es)  =  e  d’  =  esEe’  =  d  cE  d’ 

(ps  -PT)  ■  gh'  =  h  sD  h'  =  gcDg' 

These  relations  show  that  the  presence  of  piezoelectricity  modifies 
the  elastic  and  dielectric  constants,  that  the  mechanical  (elastic)  conditions 
must  be  specified  for  dielectric  quantities,  and  that  electric  conditions  must 
be  specified  for  elastic  quantities.  To  obtain  an  order  of  magnitude  estimate 
of  the  size  of  the  effect,  we  may  neglect  the  matrix  nature  of  the  above 
equations  and  treat  them  as  scalars.  Then,  using  the  symbol  A  to  represent 
differences,  we  have,  from  the  first  difference  relation: 

(Ac/c)  =  e  ps  e/c  =  e2  /(e  c)  =  k2, 

and  similarly  for  the  other  relations.  Thus  the  fractional  differences  are 
determined  by  dimensionless  combinations  of  elastic,  piezoelectric,  and 
dielectric  constants.  We  call  each  combination  the  square  of  a  quantity,  k, 
known  as  a  piezoelectric  coupling  factor.  See  Appendix  4. 


Appendix  4.  Electromechanical  Coupling  Factors 

Electromechanical  coupling  factors  k  (0  <  k  <1)  are  dimensionless 
measures  of  efficacy  of  piezoelectric  transduction,  and  are  far  more 
important  than  the  piezoelectric  constants  taken  by  themselves.  They 
appear  in  considerations  of  bandwidth  and  insertion  loss  in  transducers 
and  signal  processing  devices,  in  location  and  spacing  of  critical 
frequencies  of  resonators,  and  in  electrical/mechanical  energy  conversion 
efficiency  in  actuators.  For  high  frequency  plate  resonators,  coupling 
factors  have  the  generic  form  k  =  e  /  V(ec);  for  low  frequency  bar  and  rod 
vibrators,  the  form  becomes  k  =  d  /  V(es);  alternative  forms  using  the  “h” 
and  “g”  piezo  parameters  are  used  as  well.  These  quantities  are  also  called 
piezocoupling  factors. 

The  passage  of  an  acoustic  wave  in  a  piezoelectric  medium  produces 
a  self-generated  electric  field  that  alters  the  effective  elastic  coefficient 
determining  the  wave  velocity.  The  effective  elastic  coefficient  is  called  the 
stiffened  elastic  constant  c;  it  is  computed  from  the  cE,  e,  and  es  tensors; 
treating  these  as  scalars  for  the  sake  of  simplicity,  c  takes  the  form  c  =  cE+ 
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[•] 


c  crys 

ferroelectric;  the  material  originally  is  macroscopically  isotropic  at  all 
emperatures  because  of  the  random  orientations  of  the  crystallites.  When 
pooled  through  the  Curie  point  in  the  presence  of  a  static  electric  field  (a 
process  called  “poling”),  the  domains  are  preferentially  oriented  and  frozen 
nto  a  configuration  where  the  body  is  macroscopically  piezoelectric.  It 

. !iotropy  is 

direction.  The  result  is  cylindrical  polar  symmetry 
transverse  isotropy),  represented  by  the  symbol  oomm;  as  far  as  the  elastic, 
piezoelectric,  and  dielectric  matrices  are  concerned,  the  schemes  of 
coefficients  are  identical  to  those  for  the  class  6mm.  In  fact,  for  all  tensors 
ip  to  and  including  rank  five,  the  schemes  of  coefficients  are  identical  to 


THIN  FILM  RESONATOR  OVERMODED  RESONATOR 

ELECTRODES  PIEZOELECTRIC  ELECTRODES 
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a) .  A  thin-film  fundamental  mode  resonator  supported  at  its  periphery  by  a  substrate. 

b) .  An  overmoded  resonator  consisting  of  a  thin-film  piezoelectric  transducer  fabricated  on 
low-loss  substrate  having  reflecting  surfaces. 
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d).  A  stacked  crystal  filter  having  two  layers  of  piezoelectric  films.  Two  stacked  single  sections 
are  shown  connected  in  series  so  as  to  bring  the  output  electrode  to  the  surface.  (After  Ref.  [6]). 
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Figure  2.  Resonator  equivalent  networks,  a).  The  traditional  Butterworth-Van  Dyke  (BVD)  circuit 
for  lateral  excitation  (LE)  of  a  plate  resonator. 
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piaie  resonator.  s  ne  negative  capacitor  arises  Iron 
electric  field  of  the  acoustic  wave  and  the  imposed 
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Hyperbola  :  r  =  cjc} 


Figure  4.  Normalized  input  capacitance  of  the  BVD  circuit  of  Fig.  3  a; 
function  of  frequency,  normalized  to  the  resonance  frequency.  The 
hyperbolic  region  determines  the  capacitance  ratio,  r.  The  quality 
factor,  Q,  may  be  found  from  either  the  extrema  of  the  curve,  or  from 
the  slope  maximum  in  the  region  of  anomalous  dispersion. 
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rcuix  oi  rsg0  o  iu  muud  a  comf 

tic  transmission  line  represents  the  n 
5  connected  to  add  to  the  surface  trac 
mechanical  port  of  the  driving  resonator.  This  circuit  can  be  used  to 
model  thin  piezo-films  applied  on  a  substrate,  or  to  determine  the 
properties  of  piezoelectrically  inert  materials.  Z1s  Z2  are  the 
mechanical  impedances  of  the  media;  0i,  02  are  the  corresponding 
phase  “distances”  of  each  structure  (wavenumber  times  length). 
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Figure  11.  Static  and  motional  time  constants,  representing  loss  mechanisms,  and  appearance  in 
BVD  circuits.  Internal  friction  in  the  resonator  material  is  represented  by  R1  in  series  with  the 
motional  arm,  and  calculated  from  the  material  viscosity.  Losses  arising  from  ohmic 
conductivity  and  dielectric  loss  are  subsumed  into  an  effective  conductivity,  and  appear  as  a 
shunt  resistor  Rn  =  1/Gn. 
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Figure  12.  Input  admittance  of  an  open-circuited  transmission  line.  Its  distributed  nature  is 
reflected  in  the  transcendental  tangent  function  describing  its  behavior.  The  tangent  function 
expansion  in  partial  fractions  about  its  singularities  is  valid  for  all  nonsingular  points,  and  is 
used  to  arrive  at  a  description  of  the  transmission  line  with  loss  in  terms  of  a  sequence  of  RLC 
lumped  elements,  one  for  each  singularity,  i.e.,  harmonic  resonance. 


plate  resonator  with  electrodes  of  such  a  thickness  that  wave 
transmission  effects  cannot  be  neglected.  St  is  usually  the  case  that 
the  electrode  can  be  reduced  to  a  lumped  inductor  representing  a 
purely  inertial  effect. 


mm 


Figure  14.  High  frequency  ceramic  resonator  equivalent  circuit  for  thickness  excitation.  The 
series  RLC  arm  represents  the  plate  resonance  at  the  Mth  harmonic.  CR  models  the  effect  of  the 
nonresonant  modes,  and  is  particularly  important  for  high  coupling  materials.  Electrode  inertia  is 
represented  by  Le,  while  G0  models  dielectric  and  ohmic  losses.  C0  is  the  limiting  capacitance  of 
the  plate  at  very  high  frequencies. 


Q,  vs  capacitance  ration,  r  =  C0/C.,,  for  various  classes 
ty  F  in  the  figure  is  the  figure  of  merit,  F  =  Q/r.  When  th 
3  ceases  to  intersect  the  Gp  -  f  plane,  and  the  resonanc 
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susceptance  vs  conductance 


Figure  16.  Network  functions  for  the  circuit  of  Fig.  14  with  the  values  employed  in  the  numerical 
example,  and  aluminum  electrodes,  unless  otherwise  noted,  b).  Input  capacitance  vs  frequency 


f  (hertz) 

Figure  16.  Network  functions  for  the  circuit  of  Fig.  14  with  the  values  employed  in  the  numerical 
example,  and  aluminum  electrodes,  unless  otherwise  noted,  c).  Phase  angle  at  the  input  vs 
frequency 
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Figure  16.  Network  functions  for  the  circuit  of  Fig.  14  with  the  values  employed  in  the  numerical 
example,  and  aluminum  electrodes,  unless  otherwise  noted,  d).  Absolute  admittance  versus 
frequency  for  the  resonator  with  electrodes  of  aluminum,  chromium,  and  gold.  The  frequency 
shifts  are  proportional  to  the  mass  densities,  and  must  be  accounted  for  in  any  accurate 
description  of  the  resonator. 


characterizing  modern  ceramic  resonators. 
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Figure  17.  Netw« 
advanced  lumpe 
has  been  constri 
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TABLE  I. 


TABLE  II. 


TABLE  III. 


DIELECTRIC  PERMITTIVITY  [s]  MATRIX  FOR  POINT 
GROUP  6mm  AND  TRANSVERSE  ISOTROPY. 

E-fi  0  0 

o  e„  o 

0  0  833 


PIEZOELECTRIC  [e]  MATRIX  FOR  POINT  GROUP  6mm 
AND  TRANSVERSE  ISOTROPY. 

0  0  0  0  e15  0 

0  0  0  @15  0  0 

®31  ®31  ®33  0  0  0 


ELASTIC  [c]  MATRIX  FOR  POINT  GROUP  6mm 
AND  TRANSVERSE  ISOTROPY. 


C11 

@12 

C13 

Ci2 

C11 

C13 

c13 

C13 

C33 

0 

0 

0 

0 

0 

0 

0 

0 

0 

C66  =  (C11  "  C12)/2 


0  0  0 

0  0  0 

0  0  0 

C44  0  0 

0  c44  0 

0  0  c66 
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TABLE  IV.  VALUES  OF  PARAMETERS  OF  COMPLETE  CIRCUIT, 
ALUMINUM  ELECTRODES. 


Parameter 


|Yin|=  max 
Gp  =  max 
Bp  =  max 
Bp  =  infl 

dBp/dco  at  Bp  =  infl 


dCp/dco  at  Cp  -  infl 
Cp  =  min 

Cp  =  0;  dCp/dco  <  0 
Cp  =  0;  dCp/dco  >  0 


Value 

|Yin|  =  46.385304382  S 
Gmx  =46.383632396  5 
Bmx  =23.470063282  5 
Binf)  =  0.257642171  S 
=  -  39.34874044 
Bmn  =-22.913084113  5 
Cmx  =  17.34048451  nF 
Cinf,  =  0.2201134187  nF 
=  -  29.04709396 
Cmn  =- 16.89953315  nF 


Frequency 

f  =  215.6001009MHz 
f=  215.6012272MHz 
f=  215.4137051MHz 
f=  215.6013114MHz 
^s/ohm 

f  =  215.7889172MHz 
f=  215.4135398MHz 
f=  215.6011475MHz 
fF  -  second 
f  =  215.7887560MHz 
fR  =  215.6023537MHz 
fA  =  244.9161150MHz 


TABLE  V.  VALUES  OF  PARAMETERS  OF  CIRCUIT  WITH  Le  =  0, 

CORRESPONDING  TO  AN  ABSENCE  OF  MASS  LOADING. 


Parameter 

|Yinj=  max 
Gp  =  max 
Bp  =  max 
Bp  =  infl 

dBp/dco  at  Bp  =  infl 
Bp  =  min 
Cp  =  max 
Cp  =  infl 

dCp/dco  at  Cp  =  infl 
Cp  =  min 

Cp  =  0;  d Cp/doo  <  0 
Cp  =  0;  dCp/dco  >  0 


Value 

|Yjn|  =  45.695069046  S 
Gmx  =45.693288121  5 
Bmx  =23.131662750  5 
Binfl  =  0.264757835  S 
=  -37.063308726 
Bmn  =-22.561117558  5 
Cmx  =  16.71462358  nF 
Cinfl  =  0.2206966827  nF 
=  -26.75779408 
Cmn  =-16.27336196  nF 


Frequency 

f  =  220.4523222  MHz 
f  =  220.4535472  MHz 
f  =  220.2574248  MHz 
f=  220.4536353  MHz 
^is/ohm 

f  =  220.6498443  MHz 
f  =  220.2572478  MHz 
f  =  220.4534596  MHz 
fF  -  second 
f  =  220.6496720  MHz 
fR  =  220.4547724  MHz 
fA  =  249.9119070  MHz 
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TABLE  VI.  RESONANCE  AND  ANTIRESONANCE  FREQUENCIES  FOR 
VARIOUS  ELECTRODE  MATERIALS. 


Electrode  material  fR  (Mhz) 


fAiMhz)  afA^fRHiH%) 


massless  electrode 
Aluminum 
Chromium 
Gold 


220.4547724 

215.6023537 

208.0472644 

190.6250603 


249.9119070 

244.9161150 

237.0443003 

218.4866605 


11.78285384 

11.72550452 

11.59881436 

11.14464008 


TABLE  VII.  VALUES  OF  PARAMETERS  OF  CIRCUIT  WITHOUT  -C0  and 
-G0.  RESULTS  APPROPRIATE  TO  LATERAL  EXCITATION; 
ALUMINUM  MASS  LOADING. 


Parameter 

|Y|nl=  max 
Gp  =  max 
Bp  =  max 
Bp  =  infl 

dBp/dco  at  Bp  =  infl 
Bp  =  min 
Cp  =  max 
Cp  =  infl 

dCp/d©  at  Cp  =  infl 
Cp  =  min 

Cp  =  0;  dCp/doo  <  0 
Cp  =  0;  dCp/d©  >  0 


Value 

|Yin|  =  42.069823250  S 
Gmx  =  42.067454235  S 
Bmx  =  21.349182123  S 
BinfI  =  0.298882071  S 
=  -35.54129966 
Bmn  =  -  20.717786073  S 
Cmx  =  13.88404117  nF 
Cinfl  =  0.21533716  nF 
=  -  23.09600438 
Cmn  =-13.45271446  nF 


Frequency 

f=  244.9159862  MHz 
f=  244.9173998  MHz 
f=  244.7290990  MHz 
f=  244.9174751  MHz 
ns/ohm 

f=  245.1058507  MHz 
f=  244.7289519  MHz 
f=  244.9173296  MHz 
fF  -  second 
f=  245.1057080  MHz 
fR  =  244.9188136  MHz 
fA  =  268.8881905  MHz 


57 


TABLE  VMS.  VALUES  OF  PARAMETERS  OF  CIRCUIT  WITH  CR  and  Le  =  0. 


Parameter 

|Yinj~  max 
Gp  -  max 
Bp  -  max 
Bp  -  inffl 

dBp/dco  at  Bp  =  infl 
Bp  =  min 
Cp  =  max 
Cp  =  Inf! 

dCp/d©  at  Cp  =  infl 
Cp  =  min 

Cp  =  0;  dCp/dco  <  0 
Cp  =  0;  dCp/dco  >  0 


Value 

|Yin|  =  41.302065813  S 
Gmx  =41.300259872  3 
Bmx  =  20.923000271  S 
Binf|  =  0.254698334  S 
=  -  33.49994878 
Bmn  =-20.376776814  8 
Cmx  =  15.01346206  nF 
Cinfl  =  0.2089434853  nF 
=  -24.01707857 
Cmn  =-14.59569168  nF 


Frequency 

f=  221 .9957814  MHz 
f=  221 .9970788  MHz 
f=  221.8009557  MHz 
f=  221.9971664  MHz 
fis/ohm 

f=  222.1933755  MHz 
f=  221.8007800  MHz 
f=  221.9969918  MHz 
fF  -  second 
f=  222.1932044  MHz 
fR  =  221.9983766  MHz 
fA  =  249.9118350  MHz 


TABLE  IX.  VALUES  OF  PARAMETERS  OF  CIRCUIT  WITH  CR  =  0. 


Parameter 

|Yin|=  max 
Gp  =  max 
Bp  =  max 
Bp  =  infl 

dBp/dco  at  Bp  =  infl 
Bp  =  min 
Cp  =  max 
Cp  =  infl 

dCp/dco  atCp  =  infl 
Cp  =  min 

Cp  =  0;  dCp/d(o  <  0 
Cp  =  0;  dCp/dco  >  0 


Value 

|YJ  =  41.302005575  3 
Gmx  =41.300270484  8 
Bmx  =20.917603958  S 
Binf,  =  0.249734598  S 
=  -  34.86793248 
Bmn  =-20.382202679  S 
Cmx  =  15.31271378  nF 
Cjnf|  =  0.20894304  nF 
=  -25.5031160216 
Cmn  =  - 14.89494338  nF 


Frequency 

f=  217.5974847  MHz 
f=  217.5987066  MHz 
f=  217.4102764  MHz 
f=  217.5987887  MHz 
|is/ohm 

f=  217.7873002  MHz 
f=  217.4101108  MHz 
f=  217.5986246  MHz 
fF  -  second 
f=  217.7871392  MHz 
fR  =  217.5999287  MHz 
fA  =  244.9604595  MHz 


TABLE  X.  FREQUENCIES  FOR  WHICH  |Yin|  =  MAXIMUM. 


Circuit  condition 

Complete  circuit 
-C0  =  -G0  =  0 
*CR  =  Le  =  0 
Le  =  0 
CR  =  0 


Frequency  (MHz) 

215.6001009 

244.9159862 

221.9957814 

220.4523222 

217.5974847 


|YJ  .(siemens) 

46.385304382 

42.069823250 

41.302065813 

45.695069046 

41.302005575 
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TABLE  XI.  FREQUENCIES  FOR  WHICH  Gp  =  MAXIMUM. 


Circuit  condition 

Complete  circuit 
-Cq  “  -Go  =  0 
Cr  =  Le  =  0 
Le  =  0 
Cr  =  0 


Frequency  (MHz) 

215.6012272 

244.9173998 

221.9970788 

220.4535472 

217.5987066 


Gp  (siemensl 

46.383632396 

42.067454235 

41.300259872 

45.693288121 

41.300270484 


TABLE  XII.  FREQUENCIES  FOR  WHICH  Bp  =  MAXIMUM. 


Circuit  condition 

Complete  circuit 
-Cq  =  -Gq  =  0 
Cr  =  Le  =  0 
Le  =  0 
Cr  =  0 


Frequency  (MHz) 

215.4137051 

244.7290990 

221.8009557 

220.2574248 

217.4102764 


Bp  (siemens) 

23.470063282 

21.349182123 

20.923000271 

23.131662750 

20.917603958 


TABLE  XIII.  FREQUENCIES  FOR  WHICH  Bp  =  INFLECTION  POINT. 


Circuit  condition 
Complete  circuit 

-Cq  =  -Gq  =  0 
Cr  =  Le  =  0 
Le  =  0 
Cr  =  0 


Frequency  (MHz) 

215.6013114 

244.9174751 

221.9971664 

220.4536353 

217.5987887 


Bp  (siemens) 

0.257642171 

0.298882071 

0.254698334 

0.264757835 

0.249734598 


TABLE  XIV.  FREQUENCIES  FOR  WHICH  Bp  =  MINIMUM. 


Circuit  condition 

Complete  circuit 
-Cq  =  -Gq  =  0 
Cr  =  Le  =  0 
Le  =  0 
CR  =  0 


Frequency  (MHz) 

215.7889172 

245.1058507 

222.1933755 

220.6498443 

217.7873002 


Bp  (siemens) 

-22.913084113 

-20.717786073 

-20.376776814 

-22.561117558 

-20.382202679 
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TABLE  XV.  FREQUENCIES  FOR  WHICH  Cp  =  MAXIMUM. 


Circuit  condition  Frequency  (MHz)  CB  /  Cs 


Complete  circuit 
-Cq  =  -Go  =  0 
CR  =  Le  =  0 
Les® 

CR  =  0 


215.4135398 

244.7289519 

221.8007800 

220.2572478 

217.4101108 


88.56222937 

70.90930117 

76.67753861 

85.36579969 

78.20589265 


TABLE  XVI.  FREQUENCIES  FOR  WHICH  Cp  =  INFLECTION  POINT. 


C5reuit  condition 


Frequency  (MHz)  CB  /  C£ 


Complete  circuit 
-Cq  =  -Gq  “  0 
Cr  =  Le  =  0 
Le  -  0 
CR  =  0 


215.6011475 

244.9173296 

221.9969918 

220.4534596 

217.5986246 


1.124174719 

1.099781205 

1.067127068 

1.127153626 

1.067124821 


TABLE  XVII.  FREQUENCIES  FOR  WHICH  Cp  =  MINIMUM. 


Circuit  condition 

Complete  circuit 
-C0  =  -G0  =  0 
Cr  =  Le  -  0 
Le  -  0 

Cr  “  0 


Frequency  (MHz) 

215.7887560 

245.1057080 

222.1932044 

220.6496720 

217.7871392 


CglC o 

-86.31017952 

-68.70640684 

-74.54387988 

-83.11216527 

-76.07223381 


TABLE  XVISS.  RESONANCE  AND  ANTSRESONANCE  FREQUENCIES, 
ALUMINUM  ELECTRODES. 


Circuit  condition 

Complete  circuit 
-C0  =  -G0  =  0 
Cr  =  Le  =  0 
Le  ”  0 
Cr  =  0 


!r.  (Mhz) 

215.6023537 

244.9188136 

221.9983766 

220.4547724 

217.5999287 


hMm 

244.9161150 

268.8881905 

249.9118350 

249.9119070 

244.9604595 


11.72550452 

9.58775076 

11.16538336 

11.78285384 

10.94421232 
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[Tables  for  Al,  Cr,  and  Au  have,  for  the  two  cases  each  where  Le  =  0,  the 
same  values  for  fR  and  for  fA  ,  as  they  must.] 

TABLE  XIX.  RESONANCE  AND  ANTIRESONANCE  FREQUENCIES, 
CHROMIUM  ELECTRODES. 


Circuit  condition 

Complete  circuit 
-C0  =  -Go  =  0 
Cr  =  L.  =  0 
L.  =  0 
CR  =  0 


fg  (Mhz) 

208.0472644 

237.0466749 

221.9983766 

220.4547724 

210.8100702 


!a  (Mhz) 

237.0443003 

260.9394918 

249.9118350 

249.9119070 

237.3170149 


KfA^fgHfiLm 

11.59881436 

9.55712676 

11.16538336 

11.78285384 

10.60277788 


TABLE  XX.  RESONANCE  AND  ANTIRESONANCE  FREQUENCIES,  GOLD 
ELECTRODES. 


Circuit  condition 

Complete  circuit 
-C0  —  -Gg  =  0 
Cr  =  Le  =  0 
Le  =  0 
Cr  “  0 


fg(Mhz) 

190.6250603 

218.4883918 

221.9983766 

220.4547724 

195.2873702 


-MMh.z) 

218.4866605 

241.8268696 

249.9118350 

249.9119070 

219.8428424 


KfA^fgHfil m 

11.14464008 

9.33539112 

11.16538336 

11.78285384 

9.82218888 


TABLE  XXI.  SUSCEPTANCE  AND  CAPACITANCE  SLOPES  AT  INFLECTION 
POINTS;  ALUMINUM  ELECTRODES. 


Circuit  condition 


Complete  circuit 
-C0  =  -G0  =  0 

Cr  =  Le  =  0 
Le  =  0 
Cr  =  0 


dBg/dco 

(10*6  second/ohm) 
-39.34874044 
-35.54129966 
-33.49994878 
-37.06330873 
-34.86793248 


mjCo)Mco 

(10-4  second) 
-1.483508373 
-1.179571215 
-1.226612797 
-1.366588053 
-1.302508479 
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